NPS ARCHIVE 



1969 
HEAD, H. 


AUTOMATIC SHIELD CONTROL FOR A 
MODERATE TEMPERATURE ADIABATIC 
CALORIMETER 




by 

Herbert Weldon Head 



V 



SSBS5SS5S. c-oa. 

MONTEREY. CA 93943-5101 



United States 
Naval Postgraduate School 




THESIS 



AUTOMATIC SHIELD CONTROL FOR A MODERATE 
TEMPERATURE ADIABATIC CALORIMETER 

by 

Herbert Weldon Head 



June 1969 



Tki& document kaA been approved fan public ne- 
lea&e and 6 ale; tt& dlbtnibutlon l& unlimited. 




'^s 



s. 








Automatic ShieLd Control for a Moderate 



Temperature Adiabatic Calorimeter 



by 



Herbert Weldon Head 
Captain, United States Army 
B.S., Worcester Polytechnic Institute, 1964 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



from the 

NAVAL POSTGRADUATE SCHOOL 
June 1969 






tot « f\V ' 

I Un 
V\C p • 









ABSTRACT 



An automatic, electronic controller was designed for use with 
adiabatic calorimeters. The controller employed reset action, and 
its behavior at different system gains was analyzed using the 
analytic techniques of feedback control theory. In accord with the 
analysis satisfactory performance was obtained with a calorimeter 
characterized by large thermal lags, but marginal performance was 
obtained with another calorimeter characterized by a rapid thermal 
response. A second controller, which employs reset and proportional 
action gave satisfactory results during preliminary testing for both 
calorimeters . 
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I. INTRODUCTION 



When thermodynamic measurements are made in an adiabatic calori- 
meter, the heat exchange between the sample and its enclosure should 
be zero. To accomplish this goal, it is necessary to observe and con- 
trol the temperature difference between the sample and its enclosure, 
or shield. 

A calorimetric experiment begins with the sample at a known 
temperature and in thermal equilibrium with its surroundings, the 
shield. A measured amount of energy is then added to the sample. The 
shield control system (controller) must act to re-establish thermal 
equilibrium between the sample and the shield within a reasonable 
length of time and in such a way that little, or no, heat is lost from 
the sample. After equilibrium is established, the temperature of the 
sample is again recorded. During equilibirum periods, the controller 
must supply enough heat to the shield to offset losses by conduction 
through its supporting structure and by radiation to the calorimeter 
jacket, which is maintained at some constant lower temperature by its 
contact with a thermal reservoir or constant temperature bath. During 
each period when the sample is heated, the shield requires additional 
heat in order to increase its temperature at a rate to match that of 
the sample. 

In a typical system, a differential thermocouple is utilized to 
sense the sample- shie Id temperature error, thus providing information 
by which automatic regulation of the shield temperature is accomplished. 
The error signal is thus inherently a low level signal, which implies 
that the control system have high gain. To maximize sensitivity, it 
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is of importance to avoid electromechanical components and the inherent 
deadband which they induce into the control system and to maximize the 
gain, consistent with stability, with which the error signal is intro- 
duced into the control system.^ The most convenient method of supply- 
ing and controlling the shield heat is electrical, such as by electroni- 
cally controlling a current to a resistance heater on the shield. 

Automatic thermal control devices fall into two categories: 

(1) intermittently acting devices, such as a thermostat, and (2) con- 
tinuously acting, or modulation, devices incorporated into a closed 
loop feedback system. The modulating type control action is preferred 
over the on-off action since the feedback path continually senses the 
effect of the manipulated variable, the shield current, on the equili- 
brium process. Electronic control action is commonly categorized into 
three processes, proportional, rate and reset, any combination of 

which may be incorporated into the design of a modulated shield heater 

2 

current, control system. Proportional control produces a heater cur- 
rent directly related to the magnitude of the input signal from the 
thermocouple. Automatic reset control is an integrating action which 
provides a floating correction proportional to the time duration as 
well as the magnitude of the error from the equilibrium set point. Rate 
control applies a correction proportional to the rate of change of the 
error signal and represents a differentiating action. 

^Davis , Elwood R. and Moiles, E. D., "Electronic Systems for Temper- 
ature Measurement and Control", Temperature - Its Measurement and Con - 
trol in Science and Industry , v. Ill, p.. 772-773, C. Herzfeld, ed. New 
York, Reinhold Publishing Corp., (1962). 

2 Ibid . , p. 771-772. 
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This paper describes the design and operation of an ail electronic, 
dc shield temperature regulator operating in the reset mode. The regu- 
lator operates as one element of a closed loop feedback system whose 
dynamic behavior is investigated using the analytical techniques of 
feedback control theory. The regulator is to be used in an adiabatic 
calorimeter employed in taking specific heat and electrocaloric measure- 
ments on ferroelectric crystals in the neighborhood of the ferro- 
electric transition. 
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II . CONTROLLER DESIGN AND PERFORMANCE 



Since there is interaction between the various parts of the system, 
such as the shield, controller and the sample, the system must be 
treated as a whole in order to design a satisfactory controller. A 
simple, yet effective, tool for visualizing the behavior of a system 
is the block diagram. The blocks represent components, or component 
assemblies, whose physical behavior may be described mathematically by 
algebraic functions under static conditions and by integrodi f ferentia 1 
equations under dynamic conditions. Figure 1 displays the pertinent 
variables and the interconnecting information streams of a closed loop 
temperature control system. Any difference between the set point 
signal and the output of the thermocouple is presented to the control 
system as an error signal. The output of the control system regulates 
the shield heater so as to correct for disturbances. The effect of 
this control action is sampled by the thermocouple sensor and fed-back 
for comparison with the set point, thus completing the feedback loop. 

Since an equilibrium state for the system requires no control action 
other than the maintenance of a steady state value of current to the 
shield heater, which corresponds to a null or zero signal from the 
thermocouple, the two- input system may be reduced to a simpler one- 
input system, and rearranged as shown in Fig. 2. For proper control 
action, the error signal must be presented to the controller with the 
proper sign, as indicated in the diagram, so that the loop is repre- 
sentative of a negative feedback system. For example, if the 
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FIGURE 1. General Closed Loop Control System 
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FIGURE 2. Single Input Closed Loop Control System with 
Laplace Notation 
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disturbance was such that the sample has just bden heated, the polarity 
of the thermocouple signal must activate the controller to increase the 
heating current supplied to the shield heater so that the shield temper- 
ature may increase to match that of the sample. 

For the controller discussed in this thesis, reset action was 
chosen as the method of control, instead of utilizing a combination of 
two or more types of control action, such as proportional plus rate or 
the more conventional use of proportional action by itself. An inte- 
grating type controller was chosen since a practical integrating device 
could be constructed simply and economically using solid state devices, 
and because the thermal lags associated with the typical adiabatic 
shield allows stability in the feedback loop to be achieved by a lag 
control network. The control action is accomplished by storing charge 
on a capacitor and modifying the charge strored by means of the error 
signal. The voltage appearing at the integrator output drives a cur- 
rent amplifier with the resistance heater on the shield for a load. 

Several factors relative to the interaction between the controller 
and the calorimeter shield must be considered In order to obtain 
proper reset action. While a large integrator time constant permits 
distrubances of a long duration to be integrated without saturation of 
the integrator, a large time constant may inhibit control action, re- 
sulting in overdamping. A decrease of the time constant will result 
in increased hunting, and may also lead to loss of control. Preampli- 
fication of the error signal increases the sensitivity of the controller, 



Heating of the sample and the measurement of Its temperature are 
accomplished by separate, independent circuity which is external to 
the thermal equilibrium feedback loop. 
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but too large a value of preamplification may cause increased hunting 
if the heater current is changed too quickly and too much. If thermal 
lags associated with the diffusion of heat through the shield are 
large enough, additional damping may not be necessary in the controller 
while if they are not, damping may be added to the control system in 
the form of a compensation network. If system damping is larger than 
desired, the compensation network will take the form of a lead net- 
work to lessen the damping. If thermal lags are too large, the con- 
troller cannot be made to give a rapid response. 

Thus, a simple reset action controller, composed of a preamplifier, 
a lead or lag compensation network, an integrator, and a current ampli- 
fier, behaves in a quantitatively complicated fashion. In designing 
the controller, the techniques of automatic control theory, to which 
the block diagram has provided a brief introduction, provide some in- 
sight into a qualitative, design solution for the performance of the 
temperature controller in a system whose thermal constants are known. 
Final values of controller components may be selected after suitable 
operating experience has been obtained using the engineering estimates 
arrived at as a result of the control theory analysis. 

It is easier to manipulate the Laplace transforms of the blocks in 

Fig. 2 than it is to work in the time domain. Thus, the following 

analysis uses the rules of transform algebra found in any standard 

4 5 

text on automatic control theory. ’ The two basic assumptions in the 

4 

Clark, R. N., Intr oduction to Automatic Control Systems, Wiley, 

1957. 

"Van Valkenburg, M. E., Network Ana Lysis , Prentice-Hall, 1955. 
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analysis are (l) the system elements can be described by lumped-para- 



meters, as opposed to distributed parameters, and hence the system may 

be described by ordinary linear differential equations rather than 

partial differential equations whose transforms are more complex, and 
\ 

(2) the response functions of system elements are linear. Non-linear 
functions are linearized by assuming the function can be approximated 
by just the first-order term of its Taylor series expansion. With 
reference to Fig. 2, one can define the following quantities in trans- 
form notation: 

G(s) is the thermocouple transfer function; 

H(s) is the controller transfer function; 

E(s) is the Laplace transform of the error, e(t); 

B(s) is the Laplace transform of the primary feedback signal, 

b(t) ; 

C(s) is the Laplace transform of the actuating signal, c(t)and 

Q(s) is the Laplace transform of the input disturbance, q(t) . 

These variables are related by 

E(s) = G(s) x C(s) , 

B(s) = H(s) x E (s) , 

C(s) = Q (s) - B (s) . 

These equations may be combined to eliminate B(s) and C(s), 



E 1*1 

Q(s) 



G(s) , 

1 - G(s)H(s) 



W(s) 



where W(s) is the system transfer function. 

The thermocouple has a first order time constant response, since it 
is basically an energy storage element in series with a resistance to 
flow. The forward path transform has the following form: 



G 



c 

s + d 
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The controller transform, H, of Fig. 2 is constructed from the 
transforms of the individual controller stages - preamplifier, compen- 
sation network, integrator and buf fer/amplifier . A high gain preampli 
fier with adjustable gain, , has a small time delay, 1/e, in its 
response so that its transform is 

eK. 

c - A 

°A s+e * 

A first-order phase lead compensation network has the form 

s + a . , , 

G = — - — — with b > a 

c s + b 



An ideal operational amplifier, which is used to construct an analog 
integrator, has gain, input impedance, and bandwidth which approach 
infinity and zero output impedance. Although real operational 
am plif iers have limitations, which complicate the form of the inte- 
grator transform, currently available, solid state, dc operational 
amplifiers have characteristics which approach the ideal so that the 
complete integrator transform may be approximated by a simple function 



G t « -1/RC . 

I s 

The buffer and current amplifier stages following the integrator have 
transforms which are pure gains taking the forms and -K , respect- 
ively. If a stage does not load the previous stage, the complete 
controller transform can be written 

H = g.g ; g t k b (-k ) = 



with 



Ac 



K = K. K K /RC 
A B p 



sec 



s (s + b) 
1 



Tobey, G., "Analog Integration", 
v. 42, p 49-50, January 1969. 



Instruments and Control Systems , 
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The system transfer function may now be specified as 



W = 



c (s+b) (s+e) s 



1 pu 4 3 2 

S + (b+d+e)S + (be+bd+de)S + (bde + Kce)S + Kcea 



where, the open loop gain is given by 

rH = Kec (s+a) 

S(s+b) (s+d) (s+e) * 

and 1 + GH is defined to be the system characteristic equation. Figure 
3 is a root locus plot of the characteristic equation which shows the 
loci of the roots in the s plane as the gain, K, increases from zero 
toward infinity. For K = 0 all the roots lie in the stable, left-half 
of the s plane. As K is increased toward infinity, the poles move to- 
ward the zeros (there are 3 zeros at minus infinity) along the dipole 
branches and the pole pair near the origin enters the complex plane as 
a conjugate pair. The imaginary axis is a region of marginal stability 
and the system will exhibit "ringing" at a critical value for K which 
places poles on this axis. From Fig. 3, one sees that oscillation will 
occur sooner for large K in a system with small thermal constants, c 
and d, than in a system in which the shield- thermocouple arrangement 
has a large thermal lag. 

The system behavior is examined by applying a disturbance, Q, in 
the form of a step signal, T/s, and examining the resultant response, 



E = QW = 



Tc (s + b) (s + e) 



4 3 2 

S + (b+d+e) S + (bc+bd+de)S + (bde + kce)s + Kcea 



The general procedure is to factor the denominator into four roots, 
expand E using partial fractions, solve for the residues and perform 
the inverse transform to the time domain with the aid of a table of 
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locus 
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FIGURE 3. Loci of the Roots of the Characteristic Equation 
as a Function of K 



Laplace transforms. The solution of the quartic with particular numer- 
ical coefficients is straightforward, albeit tedious, but it is too 
difficult to attempt for the general case. However, certain approxi- 
mations may be made which reveal the nature of the solutions for two 
cases; K small, say K~ 10, and large K in the region where ringing 
occurs . 

For small K (and small s), the pole due to the preamplifier lag is 
considered 'remote' by virtue of its location in the remote part of the 
s plane. The residue at this pole, is small compared to the other 
residues and, in addition, the time constant of this pole is very short 
when compared to the other time constans so that its contribution to 
the system behavior is negligibly small. For small K, with ~ 
and Q = T/s, the error behaves like 



At the output of the termocouple, this equation has an initial value of 



The disturbance is converged to zero, and since q(o) = T, by defini- 
tion, there exists a requirement that 



For large values of K, at large s, the pole of the pole-zero dipole 
branch is nearly cancelled by the zero and the residue of this pole 
becomes negligibly small. In this case, the compensation network has 
no effect, and the system behavior is characterized primarily by the 



Burington, R. S., Handbook of Mathematical Tables and Formulas , 
pp. 7-9, Handbook Publishers, 1957. 




\ 



q (0) = K e (o) = T (2b - d)/3. 



b ~ (d + 3)/2. 
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remote pole and the conjugate pole pair. Shunting of the integrator 
input by a larger b/a ratio prevents extension of the compensation net- 
work to very high gains. No attempt was made to derive the general 
solution for e(t); however, the error does converge to zero since the 
behavior of the system is the sum of a dying exponential from the re- 
mote pole and a sine term multiplied by a dying exponential from the 
conjugate pair. This behavior is valid until K is large enough to 
move the pole pair across the imaginary axis into the region of in- 

g 

stability. An application of the Routh Criterion for stability 
shows the critical gain to be 




As K is increased beyond the critical value, the system nonlinearities 
assume a dominate role. At a high enough value of K, the system re- 
sponse will ultimately be bounded by a limit cycle operation describ- 
able by rapid shutoff. Thus, automatic shield control is accomplished 
through a combination of continuous and intermittent control as K is 
increased through several orders of magnitude. 

Using the results of the previous analysis, construction of the 
controller procedes in a straight-forward manner. For measuring and 
amplifying the output of the differential thermocouple, which is 
typically fractions of a microvolt, a Keithley model 149 millimicro- 
voltmeter was chosen. The model 149 is highly stable, has nanovolt 
sensitivity and contains built-in zero suppression to cancel spurious 

potentials generated along the thermocouple transmission line and at 
the amplifier input. The amplifier output drives the controller shown 

g 

Van Valkenburg, oja. ci t . , p. 417. 



19 



schematically in Fig. 4. As the discharge of the capacitor will 
normally be controlled by the input impedance of the buffer, it is 
then required that the leakage resistance of the capacitor exceed the 
input impedance of the buffer. This requirement can be met by a 
variety of capacitors. Since a large capacitance is necessary to pro- 
duce a reasonable time constant, a tantalum capacitor was used. The 
capacitor may not be reduced with reduction compensated by an increase 

in , as the input: current offset, drift and noise of the operational 
9 

amplifier Al are proportional to and these effects may become 

troublesome enough to require separate compensation if R^ is further 

increased. ^ The output stage can deliver 4W of continuous power to 

a load drawing 140 ma, when a heat sink is used with is a 10<fl> 

temperature sensitive device which stabilizes during sustained 

11 

operations at higher currents . The controller was mounted on a 
small etched circuit plug board. A standard 5 in x 19 in rack panel 
contains the controller, the necessary power supplies, an output cur- 
rent monitor, and a battery source for manually driving the controller, 
The input is Zener limited to a + 10.2 volts to protect Al from 
overloading. 



9 



Type uA 741, Fairchild Semiconductor, 313 Fairchild Drive, 



type ua /4L, r 
Mountain View, Cali 



fornia . 



10 . 



Weden, C. V., "Review of Operational Amplifier Principles," 
Application Bulletin - July 1966, Fairchild Instrumentation, 475 
Ellis Street, Mountain View, California 

^"Sensi tor ", manufactured by Texas Instruments Inc., Post Office 
Box 5012, Dallas, Texas 
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The controller constants may be computed from the component values 

in Fig. 4. The output stage voltage gain, Kp, is derived using an 

12 

approximate two parameter equivalent current. 



K 

P 



-bfA 



h . 
le 



+ (1 + h c )R ~ 
fe e 




h. /h 
le 



fe 



+ R 

e 



1.8 for R = 200Q 
Li 



The transistor parameters were measured with a Tektronix Model 575 



Curve Tracer, and found to be h. ^ 17. 5K and h 



175. The buffer 



le f e 

stage voltage gain expression is derived from the JFET equivalent 
13 

circuit and the gain , is computed using the transistor parameters 
given by the manufacturer's data sheet, 



HR ' 

K_ = t D i T in = 0.9 where R' = R //[R, + h. + (1 + h )R ]. 

B r , + R (M» + 1) s s L b le fe e J 



The Keithley time constant was measured by observing the output volt- 
age rise time to a step input on the 100 p,V scale with a Tektronix 
model 454 CRO and found to be 

1/e =» 2 msec 

For the compensation network used 

a = 1/R^C^ = 1 sec 3 and b = l/CR^Z/R^) - 1.2 sec 

If the thermal constants for a particular system are known, the 
idealized behavior for an assumed input may be predicted from E = QW . 
An examination of the orders of magnitude of the various terms in the 



1 2 

Mi liman, J. and Halkias, C. C., Electronic Devices and Circuits , 
pp. 340-341, McGraw-Hill, 1967. 

13 Ibid., pp. 404-405. 
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FIGURE 4. Controller Circuit Diagram 




transfer function reveals which element (s) of the controller dominates 
its behavior for a given gain. The critical gain for stability of the 
system can only be calculated if the thermal delays of the calori- 
meter are specified. 




FIGURE 5. Controller Etched Circuit Board 



/ 
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III. RESULTS 



The controller was tested with two calorimeters of very different 
thermal design. The first calorimeter, for which the most extensive 
testing was carried out, was for use near 120°K and was characterized 
by sluggish thermal response. The second, a developmental calori- 
meter, was designed for use near room temperature and was character- 
ized by very small thermal time constants. 

The lower temperature calorimeter contained a cylindrical Cu 

3 

shield, with dimensions 2.0 in o.d. x 0.016 in wall x 4 / 16 in length, 
and mass of 98.7 gms . The shield was isolated from the calorimeter 
jacket by a cyclindrical Cu furnace with dimensions, 2^/8 in i.d. x 
0.002 in wall x 4^/ 16 in length, and mass of 38 gms. Both cylinders 
screwed onto a 47 gm Cu lid which was silver brazed to a % in i.d. 0.01 
in stainless- s teel tube 8.5 cm long. The entire assembly was sealed 
within a vacuum jacket which, in turn, was immersed in a bath of liquid 
nitrogen. The heater consisted on an 80 0 manganin coil around the 
lid stem in series with a 27 Q manganin coil around the outside sur- 
face of the furnace cylinder so that most of the heat appeared at the 
lid. Current distribution between the two heater coils was externally 
controlled by a potentiometer to trim gradients and the entire heater 
assembly appeared as a 200 Q load at the controller terminals. The 
thin walled stainless-steel tube was characterized by a thermal 
resistance , 

R. , = AT/Q = ^/AK = 625 °K/w 
th 
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where Kss = 0.17 W/cm - C K. A. t T ~ 120^K radiation was relatively 



small . 
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The thermal lag of the calc rime, ter appears in the system transfer 
function of the previous chapter through G, the thermocouple transfer 
function, which is d-termin=.tJ experimentally . The. testing technique 
measures the transient response of the th-cmocouple by comparing the 
measured temperature reading to a theoretical ramp change in tempera- 
ture produced by imposing a stc.p change in Current through the shield 
heater under adiabatic conditions.* - * 3 Neglecting the IR drop at the 
thermocoupLe-shi' Id junction, the calorimeter tested had time con- 
stants given approximately bv 

1 A « 6 sec and 1/d ss 3 sec. 

These constants are for a 0.002 in d i am chrome 1-a lumal thermocouple 
bonded to the shield with GE703L varnish. When these values are in- 
serted into the previous theory for e(t), the maximum allowable gain 
for continuous control action is 



K % 10 
A 



Since this thermocouple produc-s an average of I5p.7 pc r degree of 



14 



Reese, W. and May, L. 1. , p hys. Rev . 1_62_, 512 (1967). 

^Reese, W. and May, I. I 1 ., Laboratory Notebook, .January 1966 - 
June 1967, Naval Postgraduate School, Monterey, California 93940, 

unpublished . 

^Green , S. J. and Hun* , T. W. , "Accuracy end Response of the 
Thermocouples for Surface and r luid 1 emporaturt. Measur -mtnts , " 
Temperature -Its M. ? .sur ment and Control i n S c ic-nc.- and industry , 
pp. 708-711. 
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differential, the Keith Ley should be operated on the 1 p-V scale where 

K = 10 7 for maximum sensitivity. At this value of K , full scale 
A A, 

ringing occurred and the controller entered a limit cycle. However, 
the performance of the controller on the 1 |lV scale, as measured by 
the thermal response of the sample, was acceptable as shown in Fig. 6 
and was an improvement over a proportional controller which was pre- 
viously used with this calorimeter.^ 7 

For too large a sample heating step, corresponding to coarsely 
spaced data points, the controller appLied maximum available power 
initiating an extended limit cycle observed as an overshoot on the 
sample recording. The maximum sample heating step was experimentally 
determined and depended on several factors; for example, the current 
gain of the output stage, the heat capacity of the sample and the 
proximity of the thermocouple junction to the sample heater. 

Drift of the sample temperature for the sample used to obtain Fig. 
6 at a point about 5°K below the ferroelectric anomaly was improved 
from 0.045°K/hr to 0.006°K/hr, measured over a 9 houi period, after 
a few minutes of careful adjustment of the Keithley zero suppress. 

The controller was less satisfactory for the second calorimetric 
system. The second system was characterized by a light mass A1 con- 
struction and coupling to the bath was by a short brass link with R^ = 

8.3°K/W. In addition, radiation coupling of the shield to the bath 

3 

was much more important (cyT ) as was the increased coupling of the 
sample to the shield. The constants, c and d, were not determined for 
this calorimeter but the time response was at least a factor of 10 

^ 7 Reese and May, Phys, Rev . 162 , 612 (1967). 
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*The vertical scale represents temperature increase and the horizontal 
scale shows elapsed time with 50 sec between divisions 



FIGURE. 6. Sample Trace of Temperature Versus Time Illustrating 
Equilibrium of Potassium Dihydrogen Phosphate for a 
.025°K Step Beginning at T = 118.151°K. 
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faster. A sample recording exhibited unsatisfactory oscillations 
tracable to the ringing of the thermocouple signal. The trace oscilla- 
tions were minimized, but not entirely eliminated, by reducing the inte- 
grator time constnat by a factor of 10, limiting the maximum current 
delivered to the load and reducing the gain, K . Despite these changes 
the system still performed only marginal Ly ; however, the system could 
be used to accrue data. 

To overcome this problem, a proportional control stage was added 
and a second controller was assembled which operated in both the reset 
and proportional modes. The test circuit is shown in Fig. 7 with 
approximate component values. This network has exhibited extremely 
good stability with the Cu calorimeter. Control action was so posi- 
tive that the thermocouple error never exceeded 2 p-V during sustained 
sample heating under maximum power. The performance of this controller 
in the second system was an improvement over the reset action con- 
troller; however, the sensitivity of the calorimeter still exceeded 
the ability of the controller to eliminate hunting on the sample trace 
at the required gain level of = 10^ under all conditions, although 
stable operation was obtained with K. = 3.3 x 10 . 
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FIGURE 7. Circuit Diagram of Controller Operating in the Combined Reset and 
Proportional Mode 
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